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The Stability of Convergence of Curve Evolutions in Vector Fields 

Junyan Wang* Member, IEEE and Kap Luk Chan Member, IEEE 

Abstract 

Curve evolution is often used to solve computer vision problems. If the curve evolution fails to 
converge, we would not be able to solve the targeted problem in a lifetime. This paper studies the 
theoretical aspect of the convergence of a type of general curve evolutions. We establish a theory for 



analyzing and improving the stability of the convergence of the general curve evolutions. Based on this 
theory, we ascertain that the convergence of a known curve evolution is marginal stable. We propose 
a way of modifying the original curve evolution equation to improve the stability of the convergence 
according to our theory. Numerical experiments show that the modification improves the convergence of 
' the curve evolution, which validates our theory. 

u' 

Index Terms 



Object segmentation, active contour, geometry of curve evolution, stability of convergence. 

I. Introduction 

Various computer vision tasks have been formulated and solved based on the formulation of the curve 
evolution equations, such as boundary detection JT), shape-from-shading ||2l and optical flow 0. In 
boundary detection, for example, the desired object region are supposed to be enclosed by the convergent 
curve produced by the curve evolution. The requirement of convergence is very necessary. Suppose that 
the curve evolution cannot converge, the evolving curve may eventually disappear or we would never 
obtain a solution in a lifetime. 

The curve evolutions obtained by gradient descent based functional energy minimization H] 
are globally convergent in theory [6]. Furthermore, the numerical convergence of some of those curve 
evolutions was also shown by simulation Q and theoretical analysis JH. Hence, the curve evolutions 
in this framework are convergent once formulated. Beside of the methodology of energy minimization, 
people have also tried to directly formulate the object boundaries as the solution to ad-hoc curve evolution 
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equations, such as in @ iflOl ifTTl 11121 |[T3l . However, no general technique for either determining or 
ensuring the convergence of such curve evolutions is known. 

As a main contribution of this paper, a theory for analyzing and improving the convergence of curve 
evolutions has been established. First, a condition for determining the stability of the convergence of a 
general curve evolution in a vector field F is derived. There could be various ways to improving the 
stability of a curve evolution. We also expect the improved curve evolution to behave similar to the 
original curve evolution. Then, an error bound that was extensively used in control theory is also adopted 
for measuring the similarity of curve evolutions. The similar curve evolutions are supposed to have similar 
stationary solution. Finally, we show the marginal stability of the convergence of the curve evolution of 
Equilibrium Flow proposed in [13 |. We also modify the original curve evolution equation of Equilibrium 
Flow to improve the stability of the convergence. 

The rest of the paper is organized as follows. In section JIJ we introduce background of curve evolution. 
In section |III1 we present the mathematical analysis on the convergence of curve evolution in a vector field. 
In section |IVJ we present the application of the mathematical analysis to a system of curve evolutions. 
In section [V] experimental results validates our theoretical claims. We conclude the paper in section rvTl 

II. Background 

In this section, we provide the background to facilitate the subsequent discussions in this section. A 
general formulation of curve evolution can be written as follows lfl4l lfT5l : 

= a(p, t)T(p, t) + 0(p, t)N(p, t) 
dt (1) 

= a(p, t)T(p, t) + (F, N(p, t))N(p, t), 

where a, (3 are the general form of coefficients defined along the curve, F(p, t) is the underlying vector 
field which is a vector function on the 2D (image) domain, C(p, t) denotes the curve parameterized by 
p evolving with respect to time t, T is the tangent of the curve, N is the outward normal. Throughout 
this paper, we discuss about the curve evolution under the assumption below. 

Assumption 2.1: The vector field F only depends on the location of C (p, t), i.e., F(p, t) = F(x) \g^c(p,t) ■ 
This implies that the curve evolution is only due to the static vector/velocity field F on the domain. 

A fundamental property of the curve evolution of Equation (QQ) was reported by Epstein and Gage in 
O, known as the Lemma of curve evolution. 

Lemma 2.2 (Epstein and Gage): Given the closed curve C(p,t) parameterized by arbitrary p G B at 
an artificial time t with the normal N, the tangent T of the curve, and given the geometric flow of a 
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curve evolution by 

^^l = a(p,t)f(C(p,t)) + f3(p,t)N(C(p,t)). (2) 

If P does not depend on the parametrization, meaning that is a geometric intrinsic characteristic of the 
curve, then the image of C(p, t) that satisfies Equation© is identical to the image of the family of curves 
C(p,t), parameterized by p € B, that satisfies 

d -^>=Pip 1 t)N{p,t). (3) 

This can be verified by the well-known level set methods. 

III. The mathematical analysis of the convergence of curve evolution 

The results of this section are obtained based on the Lemma of curve evolution 12.21 and the idea of 
using the zero level set of a signed distance function to represent closed curves. The proofs can be found 
in the appendix. 

Lemma 3.1: Given the curve evolution C(p,t) defined as follows, 

— = e(p,t)N(p,t), (4) 

the following holds 

d^c = d^p 2 t 1 , 

dt {k+i) dt k JN WJ- PJ 

To prove the above lemma, we need the lemma below. The proof of Lemma 13.11 is straightforward by 
following Lemma 13.21 and the chain rule of differentiation. 
Lemma 3.2: 

^(C)N(C) = 0. (6) 

The above results eventually lead us to a more practical result that is stated as follows, which can be 
proved by substituting Equation © into the Taylor series. 

Lemma 3.3: If there exists R > 0, such that V|r| < R, such that the Taylor series of C(p, t + r) at 
C(p,t) exists, then there exists a unique t), such that, 

C(p, t + r) = C(p, t) + t£( P , r)N(p, t), (7) 

where 



k-l 



k=l 
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By this lemma, we convert the curve evolution in the form of PDE to an algebraic equation to facilitate 
the analysis of curve evolutions in section [TV] 

We now analyze the convergence of the curve evolutions. The objective is to ascertain whether the 
convergence is stable. Naturally, we consider the case of unstable as a lack of convergence. Given a curve 
evolution defined by the following, 

9 -^ = {F,N)N- (9) 
The convergence of the curve evolution is stated as follows: 

lim ggM = (p,N)N = 0. (10) 

t^rr Ot 

The condition of convergence can be found by considering a small perturbation on the converged 
curve. If the perturbation finally vanishes at that converged curve, the curve evolution is stably convergent 
otherwise it is unstable. We now examine the convergence of the curve evolutions following this idea. 

The curve evolution near the converged curve C* can be represented as follows. 

C* = C(p,t*) = C(p,t)\ t=t ,_ T + r,(p,r), (11) 

where i] is a general formulation of perturbation, r = t* — t > and is small enough and monotonically 
decreasing during the evolution. The fundamental condition of convergence is shown below. 

Theorem 3.4: If J nn (C*) < 0, the curve evolution is stable,lf J n n(C*) = 0, the curve evolution is 
marginally stable \H J n n(C*) > the curve evolution is unstable, where J nn {C*) = [N T J[F]N](C*), 
J[F] is the Jacobian of vector field F. 

Accordingly, we can use J[F](C*) to examine (the stability of) the convergence of curve evolution in 
vector field. The proof, which makes use of Lemma 13.21 extensively, is deferred to the appendix due to 
its length. 

Suppose that the original curve evolution is unstable or marginally stable, a modification on the 
formulation of the curve evolution equation might be performed. By using the above theorem, it can be 
verified whether the modified curve evolution is stably converged. However, this leads to another question 
that whether the modified curve evolution maintains the function of the original one. The following 
theorem could be a fundamental answer to this question. 

Theorem 3.5: For a fixed p, such that dS j£'^ = 9C q^ + G{p,t), and 9C g^ is Lipschitz continuous 
with the Lipschitz constant L, Given that \\G(p, t)\\ < > 0, then the following holds, 

\\C(p,t) - S(p,t)\\ < \\C(p,0) - S(p,Q)\\e L(t - to) 

(12) 

L 
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The above theorem provides an upper bound of the difference between curve evolutions. If the difference 
is small, the curve evolutions will behave analogously, and we suppose the stationary solutions of them 
are also similar. We shall not include the proof here, since the proof of a more detailed and general bound 
for ordinary differential equations, related to Gronw all- Bellman inequality, is clearly stated in Chapter 3 

of da. 

IV. Improving the convergence of GeoSnakes 

A. Active contours with a tangential component 

In lfl3l . it has been noticed that the full gradients of the edge indicator, e.g. the magnitude of image 
gradient, along the boundaries can be uniformly small, i.e. both the normal and tangential components 
of the gradients of the edge indicator along the boundary are close to zero. However, the conventional 
active contours, such as 0], considers only the normal component of the gradient of the edge indicator. 
The curve evolution based only on the normal component often converges at the places where the 
tangential component is still significant. Therefore, a new active contour model termed the Geodesic 
Snakes (GeoSnakes) was proposed. The Euler-Lagrange (EL) equation of the GeoSnakes requires the full 
gradient of the edge indicator, i.e. both the normal and tangential components, to vanish at the contour 
of solution. Nevertheless, the direct curve evolution to solve the EL equation only solves the normal 
component of the EL equation, and the curve evolution can converge when the tangential component 
is still significant, which we call the Pseudo Stationary Phenomenon. We propose an auxiliary curve 
evolution equation, termed the Equilibrium flow (EF), to solve specifically for the tangential component. 
The full EL equation involving both the normal and tangential components is solved by alternating the 
direct curve evolution and the curve evolution of EE The system of the alternating curve evolution 
equations is as follows. 

d t C\ = (F, N)NN, c£(t, 0) = C%-\t, oo) (13a) 
d t Cl = (KF,N)N, C 2 fc (r,0) = C^(r, oo), (13b) 

where k = 1, 2, . . . , C?(r, 0) = C§(r, oo) = C , C(r, oo) = lim Cf (r, oo) = lim C|(r, oo) and K is 

k— >oo k—^oo 

a rotation matrix of 2 x 2, F = V<? is the gradient field of a scaler field g. Equation dl3a| ) is the original 
gradient descent flow, Equation (1 1 3bb is the proposed auxiliary flow called Equilibrium Flow. Note that 
we omit the curvature term that imposes smoothness to simplify our discussions. 
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B. Application of the theory to GeoSnakes 

In this subsection, we apply the theoretical results obtained above to the system of ( [TBI to show that the 
system of curve evolutions may not converge stably. We also show that the convergence can be improved 
by modifying the curve evolution equation in accordance with our analysis. 

For the gradient descent flow in Equation d!3at > under F = Vg, where 

J nn (c*) = (v 2 gN(c*),N(c*)), (14) 

where V 2 g is the Hessian of g. We may use the following lemma for analyzing J nn {C*). 

Lemma 4.1: Given that the curve evolution C(p,t) defined by Equation d 1 3ab near the converged 
curve C(p,t*), i.e. C(p,t* — 5t) can be expanded by Taylor series at t*, within a radius of convergence 
\6t\ < R, where t* is the time of convergence, and given that the initial curve C(p, 0) is not on minima 
of g (maxima of —g), then 

(V 2 gN,N)<0. (15) 

The proof which makes use of Theorem 13.31 can be found in the appendix. 
The above result shows that J nn < on C* , which also proves the following. 
Proposition 4.2: The curve evolution under F = Vg is stable or the worst marginally stable. 
For the EF defined by Equation (|13bl > under TZF, where 

Jnn(.C*) = (V 2 gN(C*),f(C*)), (16) 

we have the following equality at C*\ 

(Vg, N) = (Vg, ±^^} = ±|| Vg||. (17) 
I|v5|| 

This is because the gradient vanishes in the tangential direction. Hence the gradient is in the normal 
direction. We now differentiate the above to obtain the following, 

- ± ||V 5 (C( S ))|| = ±^V 2 gC s = J nn , (18) 
ds \\Vg\\ 

which proves the following. 

Proposition 4.3: If we normalize the gradient field, e.g., ||V<?|| = 1, then J nn {C*) = 0, hence, the EF 
is marginally stable. 

Let us consider the following modification of the Equilibrium Flow define by Equation (113bl) : 

^V=(K T Vg + eVg,N)N, (19) 
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where e > is chosen to be a constant. It can be shown that the above modified Equilibrium Flow have 
better property of convergence by the follows. 

Proposition 4.4: By normalizing ||V<?||, e.g. \\Vg\\ = 1, the curve evolution defined by Equation([T9l 
is stable or the worst marginally stable. 

The proof is straightforward based on the previous discussions. 

We can also show that the slightly revised Equilibrium Flow (fT9b will not lose its function by showing 
that Equation ([T9l behaves analogous to the Equilibrium Flow Equation dl3b| ) using Theorem 13.51 in 
which G = e(Vg, N)N and e can be small. Therefore, a curve evolution behaving analogously to the 
Equilibrium Flow but converging more stably is obtained. 

V. Experimental validations 

In this section, some numerical experiments are conducted to validate the theory and to "visualize" the 
case that the curves can disappear, due to failing to converge before disappear, by the system of Equations 
(1 1 3ab and (1 1 3bb - We also show that this is avoided by replacing Equation (1 1 3bb with Equation ( fT9l ). The 
images are smoothed and the vector fields are extended and smoothed by the method of Gradient Vector 
Field (GVF) [18] [19]. We set e = 0.1 in ( fT9l in all our experiments for validation of the theoretical 
claims. During the implementation of the system of curve evolution equations, each switch is performed 
if the curve length remains almost unchanged, and we run (|13ab , d 1 3bb or ([T9l , and then dl3a| ). 

We designed a synthesized pattern, shown in Figures [T] and [2] In this pattern, there are three blurred 
disks to be segmented out. The intensity in the background gradually changes. In other words, there 
are gradients everywhere in the background that may disturb the marginally stable Equilibrium Flow. 
The experiment shown in Figure Q] validated the theory that some of the curves for enclosing the disks 
finally disappear while the other disks are captured since the curve evolution is not unstable. Contrasting 
the lack of convergence of the original Equilibrium Flow, the improvements by the modified framework 
involving Equation ( fT9l ) is also shown in Figure [2] The real data shown in Figures [3] and [4] is blur and 
there is gradient inside and outside the object region. We also observe the lack of convergence of the 
original framework in Figure |3] as well as the improvements by the modified framework in Figure |4] 
Figures [T] [3] |2] and |4] from the left to right, show the initialization, the converged solution by (I13al i. the 
curve evolution by (1 1 3bt and the final result. 

VI. Discussions and conclusion 

The contribution of this paper is the theory on examining the stability of convergence of the curve 
evolution equations formulated beforehand. The presented modification of the equilibrium flow involves 



Fig. 1: The original framework for the pattern. 
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Fig. 2: The modified framework for the pattern 
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Fig. 3: The original framework for real data. 
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Fig. 4: The modified framework for real data. 



an parameter e. It is still an open problem for determining the parameters of curve evolution equation. The 
presented theory is complete without choosing the specific value of the parameter, and the experiment 
with the pre-specified e has validated the theory. 

In this paper, a theory for analyzing and improving the convergence of curve evolution was established. 
The theory was also applied to a recently proposed framework of curve evolution to ascertain the lack 
of convergence of the corresponding curve evolution and to improve its stability of convergence. More 
applications of the developed theoretical results can be expected. 

Appendix 

Proof of Lemma 1X2} Firstly, we write ^?(C) more explicitly as follows, 
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Using the relation \ ^% \\ = N and || V$|| = 1 for signed distance function $ denned by C = {x, y\$(x, y) = 
0}, we obtain the following: 

V 2 * = ^(C), (A-2) 

i.e., ^j(C) is the Hessian matrix of <3? at the position C. 
Then take derivative of ||V<3?|| = 1, we have the following: 

V||V$|| = (V 2 $)V$ = 0, (A-3) 

which completes the proof. ■ 
Proof of Lemma \3.1[ We prove by induction. 
For k=l, we have 

3C 

— = e(p,t)N(p,t). (A-4) 
Take derivative w.r.t t, we obtain the following: 

92 ° - ^ t )iY(p ) t) + e (p,t)^ = ^^t),b y lernmaE2l (A-5) 



dt 2 dt ^' ; ' dt dt 

Now we set k = n, we have 



Hence, take derivative and by lemma 

Q(n+l) C _ deW(p,t) 

dt^ n+1 ) ~ 
which completes the proof. 

Proof of Lemma \3.3[ The proof is constructive. 
We expand C(p, t + r) at C(p, t) by Taylor series as follows: 



Q(n) c de n - x (p,t) , , A „ 

-rrr^ = — n i - N(p,t). (A-6) 



JV(p,*), (A-7) 



fc=0 



By Lemma [3TTI we can rewrite the above as follows: 



X3 T k d^~ 1 h{p,t) 



C(p, t + r) = C(p, t ) + VE fcf dt (k-i) | ^(P« *) = ^ *) + T ^Pi r )^> ( A " 9 ) 
Since the Taylor expansion exists, hence unique. By setting 



fe-i 



*^EVVf' <A - 10) 

fe=i 



we complete the proof. 
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Proof of Theorem \3.4\ We expect to obtain a differential equation that defines a contraction mapping, 
we then can show that the rj approaches 0. 

Firstly, according to Taylor's theorem, the curve evolution equation can be reformulated as follows: 



~dt 



t=f 



dC_ 

lit 



t=t*-T 



d_dC 
+ T di~dt 



+ 0(r) 



(A-ll) 



t=t*-T 



This equation is the basic form for us to construct the contraction mapping. Our idea is to use r\ to 
replace the C in the equation to obtain our formulation. 
Noting that 

dt d.T 

-1. (A- 12) 



dt dr 

t = t* t = 0, and — = — 
dr dt 



Differentiating Equation (fTTI ). we can write the following: 

drj(p, t) dC 
dr ~~ ~dt 



(-I)- 



(A-13) 



t=f-T 

This relationship allow us to replace C with rj for one of the terms of Equation (IA-1 It to construct 
the contraction mapping defined by a differential equation. 

Now we focus on replacing the C with rj for the other term in Equation (IA-llt . By making use of 
the Equation ©, we have 



8_dC_ 
di~dt 

( 



t=t*-T 



d_ 

at 



(F(C),N(C))N(C) 



t=t*-T 







o / 



t=t*-T 



dr] 



[N* J[F]N](C)^(C) 



t=t*-T 



(A- 14) 



where ^ ^ = is by Lemma [ 

We now need the Taylor expansion of rj as follows. 

T](p, 0) = T](p, T) - T— + 0(t) & T)(j3, T) = T— + 0(t) 

where rj(p, 0) = 0. 

Substituting the above into Equation (IA- 14b . we obtain the following. 



(A- 15) 



~dt 



t=F-T 



d_dc 

di~dt 



+ Q(t) = -[N t J[F]N){C)t]{ P ,t) + 0(t). 



(A-16) 



t=t*-T 
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Combining the above with Equation (1A-13K we finally arrive at the following, 

= -[N T (C)J[F](C)N(C)] V (p,T) + 0(r) = - J nn (C)rj(p, r) + O(r), (A-17) 



di] 



T=f-t 



where 



Jnn(C) = [N T J[F]N](C). (A-18) 

Now we almost obtain the differential equation in terms of 77, but C = C(p,t* — r) is not convenient 
to deal with. We therefore use the following approximation: 



Jnn(C) = [N T J[F]N](C*)-T 



T? TTt dJ[F] dC, - 

N ^ & m ]N 



O(r) f» [i\f T J[F]JV](C'*), (A-19) 



where is the second order derivative of the vector field, which is relatively small comparing to 

J[F}. This treatment is valid for the case where the first order derivative does not vanish. Otherwise, it 
is categorized as the marginally stable. 

To further explain how the differential equation (IA-17) and the approximation Equation (|A- 19b can 
define a contraction mapping, and thus can be used for test of convergence, we discretize the differential 
equation and neglect the higher order terms 0{r) to obtain the following: 

rf(p) = n k -\p) -a rJnnCCV^Cp) = (1 "A t J nn {C*)) V^Hp)- (A-20) 
Then we take norm of the above to obtain the following: 

\\r] k (j>)\\ = |1 -A rJ nn (C*)\ Wrf-^W = |1 -A rJ nn {C*)\ k \\r?(p)\\. (A-21) 

Noting that r = t* — t is decreasing, then at < 0. Hence, i.i.f |1 +a T~J nn {C*)\ < 1, lim = 0, 

t— >oo 

i.e., r\ will vanish. And i.i.f J nn {C*) < 0, [1 +a i~J nn (C*)\ < 1. ■ 
Proof of Lemma \4.1\ Let us consider the function Q(t) = g(C(p,t)) for any fixed p, then by 
substituting in Equation (1 1 3ab . dQ{t)/dt = \\ (Vg, N) \\ 2 > 0. Hence, 5 is nondecreasing w.r.t t. Assuming 
C(p,t*) is the converged solution of Equation (1 1 3ab where for all St > 0, (Vg,N) = 0, we have 
that 0(C(p,t*)) > g(C(p,t* - St)). 

By Theorem 13.31 there exists £(p) such that the following relation between C(p,t* — St) and C(p,t*) 
holds: 

C(p, t* - St) = C(p, t*) + £(p)N(C(p, t*)). (A-22) 
The Taylor expansion of g at the converged solution C(p,t*), can be written as follows: 



2 



g(C(p, t* - St)) = g(C(p, t*)) + £- (V'gN, N) + 0(\£\ 2 ) I . (A-23) 



2 



t=t* 
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Since g(C(p,t*)) > g(C(p,t* - St)), we have 

(V 2 gN,N) 

which completes the proof. 



< 0, (A-24) 

t=t* 
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